Kinetic Thomas-Fermi solutions of the Gross-Pitaevskii equation 
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Approximate solutions of the Gross-Pitaevskii (GP) equation, obtained upon neglection of the 
kinetic energy, are well known as Thomas-Fermi solutions. They are characterized by the compen- 
sation of the local potential by the coUisional energy. In this article we consider exact solutions 
of the GP-equation with this property and definite values of the kinetic energy, which suggests 
the term "kinetic Thomas-Fermi" (KTF) solutions. We point out that a large class of light-shift 
potentials gives rise to KTF-solutions. As elementary examples, we consider one-dimensional and 
two-dimensional optical lattice scenarios, obtained by means of the superposition of two, three 
and four laser beams, and discuss the stability properties of the corresponding KTF-solutions. A 
general method is proposed to excite two-dimensional KTF-solutions in experiments by means of 
time-modulated light-shift potentials. 

PACS numbers: 32.80.Hd, 03.75.Hh, 03.75.Nt, ll.15.Ha, 75.10.Jm 
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I. INTRODUCTION 

A significant portion of the physics encountered in 
Bose-Einstein condensates (BECs) of ultra-cold atomic 
or molecular gases can be well described by a mean field 
approximation, which replaces the complicated many- 
body wave function by a single mean field obeying the 
famous Gross-Pitaevskii (GP) equation [ij. This equa- 
tion complements the Schrodinger equation by a non- 
linear term, which accounts for two-body interactions. 
In analogy to single particle physics, stationary states 
are determined by the time-independent GP-equation 



f + V{r)+g\ij{r)\^ V(r)=/^V(r) 



(1) 



where T = -j^A denotes the kinetic energy operator, 
fi is the chemical potential, ^(r) is the potential en- 
ergy, and 5 is a constant, which is determined by the 
binary collision cross-section in the s-wave approxima- 
tion [2|. Not only the ground state of the system is 
often described with remarkable precision by means of 
the GP-equation but also the low energyexcitation spec- 
trum, including vortices and solitons [I|, yj- Solutions 
tp of Eq. Ill) may be characterized in terms of the fol- 
lowing useful quantities: the local phase S defined via 
tjj = IV'le''^, the particle density p = \4>\'^, the current 
density j = {ih/2m){Tp\'ip* — i>*'^'ip), and the velocity 
field V = j/p =: {h/m)\/S. Recall that V • j = and 
Vxv = 0ifp7^0, whereas the latter relation indicates 
that vortex filaments can only occur at density nodes. 

In this paper we consider a class of excited stationary 
solutions of the GP-equation, which are characterized by 
a compensation of the local collisional and potential en- 
ergies up to a spatially constant term, i.e., 



V{r) = Vo-gmr)\' 



(2) 



This property is well known for approximate solutions 
derived in the so called Thomas- Fermi regime [l|, where 



the kinetic energy in the GP-equation is neglected. Here, 
however, we are interested in exact solutions V'(r) of 
Eq. lU) subject to Eq. ([2]), thus requiring that ip{r) is 
an eigenfunction of the kinetic energy operator, which 
suggests the term "kinetic Thomas- Fermi (KTF) so- 
lutions" . KTF-solutions require some degree of non- 
linearity. They arise at the boundary between the 
weak interaction regime, where the dynamics retains 
Schrodinger character, and the strong interaction regime, 
where the inherent non-linearity of the GP-equation 
dominates. A natural environment for the emergence of 
KTF-solutions are optical lattices \^W, (i-e., ultracold 
gases subjected to periodic light-shift potentials), where 
they correspond to excited states at the edge of the first 
Brillouin zone. The interest in KTF-solutions in two- and 
three-dimensional optical lattices results from the com- 
bination of their remarkable formal simplicity with the 
possibility of non-zero current density fields, which can 
acquire unconventional spatial topologies, for example, 
that of a vortex-antivortex sheet [3] . This article presents 
a study of KTF-solutions in optical lattices, concerning 
their stability and their accessibility in experiments. In 
Sec. [Uwe point out that any light-shift potential arising 
in an arbitrary monochromatic light field with spatially 
constant polarization permits KTF-solutions of Eq. ([ij . 
In Sec. Illll we present a general discussion of the stability 
of KTF-solutions, which is applied to examples of partic- 
ular interest for experiments in Sec. IIVI Here, we con- 
centrate on one-dimensional (ID) and two-dimensional 
(2D) optical lattice scenarios readily obtained in exper- 
iments by means of the superposition of two, three and 
four laser beams. For the case of the ID lattice we point 
out that the corresponding KTF-solution arises at the 
boundary between the regimes of linear Bloch bands and 
non-linear Bloch bands, characterized by additional loop 
structures [a-(lO|. In Sec. |V]we generalize the considera- 
tions of Ref. 7] showing that for any 2D KTF-solution a 
suitable time-varying light-shift potential can be found, 
in order to drive the required current density. Finally, in 



Sec. I VII this general scheme is apphed to the examples of 
Sec. El 



II. KTF-SOLUTIONS IN ARBITRARY 
LIGHT-SHIFT POTENTIALS 

KTF-solutions arise in a large class of light-shift po- 
tentials JJj] according to the following scheme: Con- 
sider an arbitrary monochromatic light field E(r, t) = 
(l/A/2)(E(r)e*'^* + E*(r)e-*"*) with E(r) = E(r) e, where 
E(r) is a complex scalar and e is a spatially constant com- 
plex polarization vector satisfying i ■ i* = 1. Maxwells 
equations imply (A -I- fc^)E(r) — with k = w/c 
and e • VE(r) = 0. The light field E(r,i) gives rise 
to a light-shift potential V{r) = -3?(a) (|E(r, t)|2) = 
— 5R(a) |E(r)p, where a denotes the complex polarizabil- 
ity of the particles, which is assumed to be scalar, and 
the triangular brackets denote the time-average over one 
oscillation cycle. We now introduce the wave function 
^{v) = ^y^{a)/gE{r), which satisfies T^ = -ErccV" with 
the single-photon recoil energy E^^c = [hkY /2m. More- 
over, Eq. ([2]) holds with Vq = and thus Eq. ([T]) is sat- 
isfied with /I — -Ercc- In brief, monochromatic light fields 
with spatially constant polarization give rise to KTF- 
solutions of the GP-equation with a particle density pro- 
portional to the time-averaged intensity. This assertion 
may be generalized to include a wider class of light fields 
admitting certain types of polarization gradients. Note 
that due to the vectorial character of the electric field, 
there is in general more than one light field yielding the 
same KTF-solution. For repulsive coUisional interaction, 
g > Q and thus 5R(a) > is required, which corresponds 
to normal dispersion, and thus a negative light shift po- 
tential, which is obtained for negative detuning of the 
light field with respect to the relevant atomic transition. 
In this case, the density maxima arise in the potential 
minima, i.e., the confining potential stabilizes the gas 
against collisional pressure. For attractive collisions and 
positive detuning, the density maxima coincide with the 
potential maxima, such that the repelling potential force 
counteracts collisional implosion of the gas. 
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HI. STABILITY ANALYSIS 

Because KTF-solutions generally describe excitations, 
a central question in regard to their physical significance 
concerns their stability. The stability of solutions -0 of 
Eq. ([T]) may be considered via the grand canonical po- 
tential K[^,] = J <Pr[r,f^e + {V - ^i)m'' + gmVA 
for 0£ = -0 -|- ex with e € R, and an arbitrary normal- 
ized wavefunction x- Use of Eq. Q and /i = -Broc yields 
^if [V'e]e=o = and 



Stability requires that 



if [0e]e=O > . 



(3) 



(4) 



Henceforth, we restrict ourselves to superpositions of 
N optical travelling waves sharing the same polariza- 
tion vector e with arbitrary amplitudes and wave vec- 
tors ki,,v G {1,...,N} with k = |ky| for all v and 
A = 27r/fc denoting the wavelength. The correspond- 
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ing KTF-solution has the form ijj = X]i/=i V'i^ ^ 
with spatially constant complex amplitudes 0^. Dif- 
ferent choices of kj^ correspond to a rich variety of pe- 
riodic and quasi-periodic light-shift potentials @, in- 
cluding triangular, hexagonal or square lattice geome- 
tries. If the lattice potential is periodic, it suffices 
that each unit cell separately satisfies Eq. Qj. In 
each unit cell we may then expand the arbitrary per- 

turbation x = v-^^Y.^ «„« ^ni,..."« O^Li e"'""'^"- 

in a Fourier-series with respect to the Bravais lat- 
tice, with V denoting the unit cell volume. Normal- 
ization requires I]„i,...,„„gz |Cni,...,nivP = 1- For the 
kinetic term in Eq. ^ we find J(Pr2x*{T — i?rcc)x 

= 2i?recE„„...,„„ez|Cni,...,n«P[(E!Il«X)^~l] with 

kjy = h^/k. This term takes values larger than zero, if the 
expansion of x comprises terms Cni,...,nN of sufficiently 
high order J2}^=i''^''^i^' refiecting the fact that high fre- 
quency perturbations possess large kinetic energies. Let 
us assume repulsive collisions, i.e., g > 0. The collisional 

term g{4'X* + 4'*xY is positive then and ^if [V'e]e=o 
becomes positive for perturbations with kinetic energies 
exceeding E'roc, which thus do not contribute to possible 
instabilities of ^^{v). If ■!/'(r) comprises Fourier terms up 
to maximally first order, we may limit the Fourier ex- 
pansion of relevant perturbations x to second order. For 
attractive collisions {g < 0) , the negative collisional term 
gii'X* +'4'*x)'^ does not permit stability. As discussed for 
specific examples below, stable KTF-solutions typically 
require gp > -Ercc, where p denotes the mean particle 
density. Expressing g = ATrh^a/m in terms of the s-wave 
scattering length a and the atomic mass m, this leads 
to a > Omin with Oinin = k'^/Snp. Inserting typical val- 
ues (fc = 27r 10^ m-\ p = 10^0 m'^) yields amin ~ 300 ao 
(ao = Bohr-radius), which is well in reach of experiments, 
if necessary by exploiting a Feshbach resonance [ij . 



IV. EXAMPLES OF KTF-SOLUTIONS IN 
OPTICAL LATTICES 

In the following, we apply the previous general con- 
siderations to three elementary examples: a ID optical 
lattice and 2D optical lattices with triangular and square 
geometries. The collisional interaction is assumed to be 
repulsive and correspondingly the required light-shift po- 
tentials are negative. 



A. ID Lattice 

We begin with a ID optical lattice composed of two 
counterpropagating travelling waves sharing the same po- 
larization, thus yielding a light-shift potential V{x) — 
—2Vsm^{kx) with the mean potential well depth V > 
0. The corresponding KTF-solution can be written as 
(l){x) = i^/2pcSm{kx), where in accordance with Eq. ([2]) 
Pc = V/g. The simplicity of this example lets us di- 
rectly see that it marks the boundary between the regime 
of weak interaction, where the description in terms of 
Bloch states yields a conventional dispersion relation, and 
the strong interaction case, where the band structure ac- 
quires additional loops at the edges and the centre of the 
first Brillouin zone, with the consequence of inherently 
non-linear hysteretic dynamics [8l-[l0|. We suspect that 
this observation generalizes to be a generic property of 
KTF-solutions in optical lattices also in 2D and 3D. To 
clarify this point we remind ourselves (following Ref. Q) 
that the Bloch states near the edge of the first Brillouin 
zone may be approximated within a two-level picture as 
4'a,Kix) = v^(cosae*'''^^+'*^^ -I- sin a 6*'"'"^+''^^) with k 
denoting the (small) deviation of the quasi-momentum 
from its value k at the zone edge (|k| ^ 1). The varia- 
tional parameter a is determined by subjecting (f>a,K{x) 
to Eq. ([T]), which leads to the system of equations 



* l'(''>'Ere 



2^rcc + (3— -2)V~2fi\ sm(2a) + y = 

Pc 



(5) 



4E,.cc K - ^y cos(2a) sin(2a) - F cos(2a) = 

Pc 



In Fig. [T] we plot the chemical potential /x(k), solution 
to Eq. ([5]) , in the vicinity of the zone edge for the three 
cases p < Pc (a), p = Pc (b) and p > pc (c). The case 
p = Pc (b), characterized by a cusp arising at the zone 
edge (k = 0), separates the regimes with (c) and without 
(a) a loop structure in the dispersion function p{k). By 
setting K = in Eq. ^ one may directly determine the 
solutions at the zone boundary. Two solutions fj,± = 
Ercc— ^{2—3-^^l)V arise for arbitrary values of p. If p > 
Pc an additional solution (the top of the loop) ^o = E^cc — 
{1~-^)V becomes possible, for which sin(2a) — h- holds. 
The corresponding wave function </'o(a;) = i^/2ps\n{kx) 
is identical with the KTF solution 0(x) at the critical 
density p = pc and the chemical potential becomes pq = 
i^rec in this case. 

Let us next inspect the stability of the KTF-solution 
4>{x). Upon making a Fourier expansion oi x — 
^-1/2 Y^^^.j^ ^, Cnc"'''^ up to Nth order on the unit 

cell [0, A] we may express ■^K['ip^]^=Q — cM(z)c as a 
bilinear form with respect to the vector c = (a, b), where 
a = 3?[(C_w, ..., Cjv)], b = 3[(C_Ar, ..., Cn)] and 
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FIG. 1: The first and second Bloch bands derived from Eq. 
([5]) are shown in (a), (b) and (c). (a) illustrates the linear 
regime (p = |pc), (b) shows the boundary between the linear 



The cusp at «: = 
/2pc sin(fca;). (c) 





M(-)(z) 



(6) 



and the non-linear regime (p = pc). 

corresponds to the KTF solution (^(a;) 

illustrates the non-linear regime (p = fpc) characterized by a 

loop structure in the lower band. For all cases V = i?rec. In 

(d) the eigenvalues of M(z) for N = 2 are plotted versus the 

mean well depth z = V /Ercc- 



with z = gpc/Eicc = V / E^^c and the symmetric matrices 

M\^/n{z) = (n^ + Z - 1) 6n.,n + jZ ((5„+2,m + '5„_2,m + 
6n^m+2 + Sn^m-2) ± jZ {dn+2,-m + On-2,-m + 0-n,m+2 + 

S-n,m-2)- Stability requires that all eigenvalues of M(z) 
exceed zero. Limiting the Fourier expansion to second 
order N — 2 yields the eigenvalues shown in Fig. [TJd). 
The lowest eigenvalue crosses zero at z = 1.22, where 
zEioc denotes the mean potential well depth. We also 
have extended the Fourier expansion up to third order 
obtaining the same result for the lowest lying branch of 
the eigenvalues in Fig. [1] thus confirming that in fact 
higher than second order terms in the expansion of x are 
irrelevant for the stability. 



B. 2D Square Lattice 

Next, we discuss a 2D optical square lattice scenario 
composed of two optical standing waves zsin(fca;) and 
zsin(fcj/)), oriented along the x- and y-axes with linear 
polarizations parallel to the z-axis. We admit an ar- 
bitrary phase lag 9 between the oscillations of the two 
standing waves. The specific case 6 = tt/2 has been dis- 
cussed in Ref. \T\ ■ Apart from the interesting additional 
physics accessible, the inclusion of other values of 9 in our 
analysis is essential for experiments, because 9 can only 
be controlled with finite precision. Light-field configura- 
tions of this type have been extensively applied in numer- 
ous previous experiments |12l |. The corresponding KTF- 
solution is tpg{x,y) = •y/p(e*^/^sin(fcx) -I- e~^^^^ sm{ky)) 
with p denoting the mean particle density. According 
to Eq. (121) the required light-shift potential is Vg{x,y) = 
— V{sm'^{kx) + sin^(fcy) -|- 2 cos(6') sin(fcx) sin(fcj/)) with 
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FIG. 2: Local particle density pg{x,y) for ^ = (a), 8 = 7r/4 
(b), and 9 = n/2 (c). Black (white) indicates low (high) 
particle density. In (d) the particle current density ]e{x,y) 
is shown for 6 / mr. The white dashed (c) and grey (d) 
rectangles shows a A/2 x A/2 sized plaquette. In all graphs 
an area corresponding to 3 x 3 plaquettes are shown. 



V = gp. Despite its formal simplicity, ^g{x^y) pos- 
sesses remarkable properties. For 9 = m: with integer 
n, ipeix, y) represents a two-dimensional array of station- 
ary solitons separated by nodal lines, where the particle 
density pe{x, y) becomes zero (black regions in Fig.ma)). 
For values oiO ^ wk the particle density nodes are points 
(cf. Fig. [2] (b), (c)) and a periodic pattern of vortical 
fluxes arises with alternating rotational sense for adja- 
cent plaquettes. This is illustrated in Fig. [5] (d), which 
shows the particle flux density for 9 7^ riTr, given by 
J6/(a;, y) = (h/m) p sm{9)W x z sin(fcx) sin(fcy). While 
jg scales with sin(0), its spatial structure is the same 
for all 9. In contrast, the spatial structure of the corre- 
sponding velocity field vg depends on 9. This is shown 
in Fig. [31 where vj is plotted for 9 = vr/lO (a), 6* = 7r/4 
(b) and 9 = tt/2 (c). If 6 tends to zero, each contour for 
some fixed value of \vg\ approaches the diagonal nodal 
line structure of the 9 — density distribution (cf. Fig. [5] 
(a)). Nevertheless, for all 9 ^ mr, one obtains 
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h 



(7) 
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FIG. 3: The square of the velocity field Vg is plotted for 
e = Tx/lO {a), 6 = 7r/4 (b) and 9 = 7^/2 (c). The same areas 
as in Fig. [2] are shown. 



determined by the healing length. We briefly discuss the 
implications for the vortex cores of the 9 = 7r/2-case, 
which have a radius on the order of i?coro ~ 1/k. Defin- 
ing the healing length as ^ = (87rap)~^/^ with a denoting 
tlie s-wave scattering length, upon use of Eq. ^ one may 
write (fc^)^ = E-ccc/gp = E-ccc/V with mean potential well 
depth V . Thus, the condition that the core size exceeds 
tlie healing length ^ < i?corc is equivalent to i?rec < V- 

The mean kinetic energy per particle Tg = 
J d^r pg^^^g/ J d'^rpg connected with the velocity field 
Vg increases from zero to its maximal value --Eroc as 
9 is tuned from zero to 7r/2. The total kinetic energy 
per particle is i?ioc independent of 9. The difference 
Eicc — Tg reflects the quantum pressure, corresponding 
to the ^-dependent degree of localization. The angular 
momentum per particle is Lg = J d^rr x jg/ J d^r pg — 
(S/tt^) h sm{9) z. The small difference between the fac- 
tor S/tt^ and unity accounts for the fact that the velocity 
fields from different vortices with opposite sense of rota- 
tion yield some cancellation. 

We may consider the stability of V'e(a^j u) following the 
same procedure as in the ID case illustrated in the con- 
text of Eq. ([6]) . Here, a Fourier expansion up to second 
order leads to a 26 x 26 matrix, which is diagonalized in 
order to obtain the (g > 0) stability diagram shown in 
Fig. m As is seen in the figure, the vortical particle fiux 
jg tends to stabilize ipglx^y). The largest stability range 
{V/E^cc ^ 3.1 ) arises for 9 — {n + ^)n corresponding to 
Fig. [5] (c) and Fig. [3] (c). For 9 = mr the flux vanishes 
and instability occurs for any potential well depth. For 
negative 5- values, 'il'g{x,y) becomes unstable for all 9. 



C. 2D Triangular Lattice 



showing that 'vg{x,y) represents a pure vortex- anti- 
vortcx lattice with vortex-filaments at positions kx, ky e 
TT Z. For g > the vortices are pinned at the poten- 
tial maxima, which correspond to the density minima of 
the lattice (black regions in Fig.[2](b),(c)), in accordance 
with results obtained for plain vortex lattices prepared in 
large scale traps and subsequently exposed to an optical 
lattice potential [ij, Hj] ■ 

For a superfluid, we expect the particle density to van- 
ish at the poles of the velocity field on a spatial scale 



As a third example, we briefly discuss a triangu- 
lar lattice composed of three travelling waves prop- 
agating within the xy-plane with linear polariza- 
tions parallel to the 2:-axis and wave-vectors k^ — 
/c{cos(27ri//3),sin(27rz^/3)},i^ £ {1,2,3}, mutually en- 
closing a 120° angle. The corresponding KTF-solution is 

'tp^{x,y) = a/p/3 (e*''^ -I- e~'''^/^cos(^/cy)). In contrast 
to the superposition of four travelling waves, there is no 
free phase parameter here. Inserting complex amplitudes 
for the three superimposed travelling waves would merely 
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Cs = \J pgjm. In the regime of phonon-like excitations 
{{hAkY /2m <^ pg), the resonance condition is approxi- 
mated by Cg = Cs- SRS can be extended to yield excita- 
tions with more complex spatial geometries. For exam- 
ple, the KTF-solution ipg^x^y) comprises the four com- 
ponents ±hk X, ±hk y in momentum space, where x, y 
denote the unit vectors in x- and y-directions. The ex- 
citation of each momentum component requires a pair 
of counterpropagating laser beams with frequencies w 
and uj + fl, i.e., in total eight beams are necessary form- 
ing a pair of crossed bichromatic standing waves. Sim- 
ilar multi-beam variants of SRS have been proposed as 
a means to excite vortex or skyrmion states in a BEC 

[Sill. 



FIG. 4: Stability diagram for the family of KTF-solutions 
tl;g{x,y) = yp (e'"/^ sin(fez) + e"'*/^ sin(fcy)) with mean 
particle density p and light-shift potential Ve(a;,j/) = 
— V (sin'^(fc2;) -|- sm^{ky) + 2cos(6') sin(fca;) sin(fcy)). 




FIG. 5: The particle density p^{x,y) (a), the flux density 
i<i{x,y) (b), and the square of the velocity field v^{x,y) (c) 
are plotted. Black (white) indicates low (high) values. The 
equilateral triangles (with A/2 side length) indicate a single 
plaquette of the lattice. 



yield a spatial shift of the resulting field ^. In Fig. [5] we 
show the particle density p<i{x,y) (a), the fiux density 
i<i{x,y) (b) and the square of the velocity field v^(a;,y) 
(c) . We encounter a situation very similar to the case of 
the square lattice a.t 9 = 7r/2 discussed in Fig. [2] (c), (d) 
and Fig. 13] (c), obtaining a vortex-anti- vortex lattice with 
vortex-filaments at the nodes of the particle density. 



V. EXCITATION OF KTF-SOLUTIONS 

In order to study KTF-solutions in experiments, a 
technique is required to excite them. A well-known 
method to produce excitations with a definite momen- 
tum in BECs is stimulated Raman scattering (SRS). Two 
intersecting laser beams with common polarizations, fre- 
quencies w and uj + fl, and fc- vectors ki, k2 are employed 
to produce a travelling intensity grating moving at a 
speed Cg = fi/Afc with Afc = |ki — k2|. This light grat- 
ing yields a corresponding moving light-shift potential, 
which can excite a velocity field, if J7 and Afc match 
with the Bogoliubov dispersion relation [1^, [l6[, which 
requires Cg = y^(hAk/2'm)^ + Cg with the sound velocity 



In the following, we discuss the use of SRS to excite 
arbitrary 2D KTF-solutions 4>{x,y) = \'ip{x,y)\e'^^'^^'y^ 
in the xy-plane, which may possess velocity fields with 
a complex geometry. Since il){x,y) solves the Hclmholtz 
equation we may consider the monochromatic light-field 
Ei(x,2/,i) = (l/V 2)[E(a:,^)e'"« + E*(x,2/)e-"-*] with 
E(a;, y) = yj g/lk{a)xl}{x, y) z and the corresponding light- 
shift potential V{x,y) = —^{a)\E{x,y)\'^ = —g\%lj(x,y)\'^. 
Consider the following bichromatic extension 
E2(x,2/,i) = (l/^/2)[E^^(a:,y,i)e*"* + E"*(a:,^/,i)e-*'^*] 
with E^^{x,y,t) = cos(ry) E(a;, y) -I- sin(?7) e*^* E*(a;,2/) 
(with some 77 e [0,7r/2]). Since Ei(x,2/,i) by definition 
solves Maxwells equation, this holds for E2{x,y,t) in 
very good approximation, if il ^ w is assumed. Note 
that in the following examples fl/u is on the order of 
10"^^. The bichromatic field £2(3;, y,i) yields the light- 
shift potential V^{x,y,t) = -^{a)\E"ix,y,t)\^ = 
V{x, y) (1 -I- sin(27?) cos {2S{x, y) - nt)). Hence, 

V^{x,y^t) is a sum of the stationary light-shift 
potential V{x, y) satisfying Eq. ([2]) and a modulation 
term V'^°^{x,y,t) = sm{2r])V {x,y) cos {2 S {x, y) - fit) 
with an experimentally adjustable modulation strength 
sin(2?7). V'^°'^ {x , y , t) provides a light-shift grat- 
ing moving according to the wave-vector field 
K{x,y) = 2\/S{x,y) directly proportional to the 
velocity field v(x,2/) of ip{x,y). Thus, temporary 
application of V"^°'^{x,y,t) should be a means to excite 
the KTF-solution ip{x,y), if hil is adjusted to match 
the energy difference between 'ip{x,y) and the ground 
state in the lattice potential V{x,y) [2l|. Since we are 
interested in solutions V'(a;,y), for which the coUisional 
energy gp is on the same order as the kinetic energy, we 
cannot directly apply the theory of Ref. 2X\ to calculate 
the excitation efficiency. An extension of Ref. [2l| to the 
regime of significant collisional interaction is a difficult 
venture, which requires further research. In case of 
a stable KTF-solution, the fact that we can provide 
the appropriate time-dependent potential, which drives 
the required current density, is nevertheless a strong 
indication, that this KTF-solution can be efficiently 
excited. 







FIG. 6: In (a) V„/2ix,y) is replotted from Fig. [gjc). White 
(black) regions indicate minima (maxima). In (b)-(j) the 
micro-rotor potential V^i2'^{x,y,t) is shown for Q,t = wtt/S 
with n G {0, 1, ..., 8}. The same area as in (a) is shown. The 
dashed white squares indicate a A/2 x A/2-sized plaquette. 



VI. EXAMPLES OF BICHROMATIC 
LIGHT-SHIFT POTENTIALS 




FIG. 7: (a) Optical set-up consisting of two nested Michelson- 
interferometers. PZT = piezo-electric transducer, M = mir- 
ror, BS — beam splitter, AOM = acousto-optic frequency 
modulator, PD = photo diode, S — servo electronics, (b) 
Beam configuration for triangular lattice. 



In the following, we illustrate the SRS-method in case 
of the example ipg{x,y) = \ipg{x,y)\e^^'>'-'^'y^ introduced 
in the second paragraph of Sec. IIVI The required 
bichromatic light field is Eg {x,y,t) = cos(r]) Ee{x,y) + 
sin(77)e'"*E;(a;,y) with Eg{x,y) = 'i{je{x,y)^/g/Wa)z. 
The corresponding bichromatic light-shift potential be- 



Vg'Hx,y,t) 



Vgix,y) + Vr''{x,y,t) with 



l/,"-d(a;,y,i) = Vgix,y)sm{2rj)cosi2Sgix,y) - ilt). The 
time-evolution of the modulation Vg^°'^ {x , y , t) is illus- 
trated in Fig. [6] for = 7r/2, indicating that it acts as 
a collection of microscopic rotors, which apply angular 
momentum with alternating sign within the plaqucttes 
of the square lattice V„/2{x,y). In the vicinity of each 
maximum of K-/2(a;,2/) (replotted in (a) from Fig. [2][c)) 
the micro-rotor term V32'^{x, y, t) provides a quadrupole 



potential rotating with alternating directions for adja- 
cent plaquettes. In (b)-(j) y^°^(a;, y, t) is shown for 
ilt = mr/S, with n E {0,1,..., 8}, thus illustrating a 
flu t = 7r/2 clockwise rotation with angular frequency 
0,11 = i^/2. Let us briefly estimate the resonance con- 
dition. The angular momentum applied to each plaque- 
tte is approximately given by m ft^ r^ , where r = A/4 
is the distance from the centre to the edge of the pla- 
quette. Excitation of vortices requires an angular mo- 
mentum of h per plaquette, i.e., mflur^ w h and thus 
hft ^ (S/tt^) 2_Ei.cc. For rubidium atoms and a con- 
venient lattice wavelength (A = 1030 nm) 51/27r — 3.5 
kHz. It has been recently pointed out that within a de- 
scription of the optical lattice with the periodic potential 
K-/2(a^,y) in terms of a Bose-Hubbard model [l^l, the 
rotor potential V^,'^'^ {x , y , t) simulates a staggered mag- 
netic field alternating for adjacent plaquettes [20] • Fi- 
nally, we note that the modulation term V^°'^{x, y, t) ob- 
tained for the triangular KTF-solution ip<i{x, y) of Fig. [S] 
(cf. Sec. IIVI) displays a similar behavior as that observed 
in Fig.|n]for the square lattice: Centered at each density 
node, microscopic rotors with quadrupolar shape apply 
angular momentum to the triangular plaquettes of Fig. [S] 
Experimentally, the generation of the bichromatic light 



field Eg (a;, y, i), appropriate for the excitation oi'ipg{x, y), 
is straightforward using the optical set-up illustrated in 
Fig. [7] (a), thus extending a method proven practica- 
ble in previous experiments [12!]. The monochromatic 
components Eg{x,y) and e*^* Eg(a;, y) are produced in 
two nested Michelson-interferometers, each with its two 
branches folded under 90° angle. Two laser beams with 
adjustable amplitudes, frequency separation and linear 
polarizations perpendicular to the drawing plane in Fig. [7] 
are used to couple both interferometers, which comprise 
piezo-electrically driven mirrors for servo control of the 
optical path length differences and thus the temporal 
phase differences between the standing waves produced 
in each interferometer branch. The requirement, that 
the spatial fields in the two interferometers are com- 
plex conjugates, is realized by setting their temporal 
phase differences to 9 and —9, respectively. This cor- 
responds to the adjustment of optical path length dif- 
ferences A^ = ± 9\/2tt. The frequency difference Vt and 
the modulation strength sin(277) are controlled by the 
acousto-optic frequency modulator shown in Fig. [7] (a). 

Providing the bichromatic light-field for excitation of 
the triangular KTF-solution ip,^{x,y) requires the su- 
perposition of six light beams with frequencies w and 
Lj + U according to the sketch in Fig. [7] (b). All beams 
share the same linear polarization perpendicular to the 
drawing plane in Fig. [7] (b). Three beams with the 
same frequency u are used to produce the light field 
E^{x,y) = y/g/W{a)iJj<,{x,y)z. The phase differences 
between these beams determine the position of the nodes 
of E<](a;, y). To ensure that the three beams at frequency 
a; + il yield the corresponding field e*^* E*{x,y), their 
phase differences need to be the same as those between 
the oj-beams. Unfortunately, because all beams with the 
same frequency travel along different paths, an interfer- 
ometric realization of this condition is not easily imple- 
mented. 



VII. CONCLUSIONS 

In summary, we have introduced a family of stationary 
solutions of the Gross-Pitaevskii equation with definite 
values of the kinetic energy, for which the local colli- 
sional energy is compensated by the potential energy. In 
view of the fact that this property is shared with the well 
known approximate Thomas-Fermi solutions, obtained 
upon neglection of kinetic energy, we have suggested the 
term " kinetic Thomas-Fermi" (KTF) solutions. Such so- 
lutions are particularly relevant in the context of optical 
lattice scenarios, where they represent excited states at 
the edge of the first Brillouin zone. Conditions for the 
stability of KTF-solutions have been discussed and a gen- 



eral method has been proposed to excite KTF-solutions 
in experiments by means of time-modulated light-shift 
potentials. We have applied our general considerations 
to a few elementary examples: a ID optical lattice, a 2D 
square lattice and a 2D triangular lattice, however, they 
should apply to more complex lattice geometries includ- 
ing quasi-periodic lattices. 
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